      SUBROUTINE SETUP
      RETURN
      END
C   HS62
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS62'/
      DATA (X(I),I=1,3)/.7D0, .2D0 ,.1D0 /
      DATA N/3/ , NH/1/ , NG/6/
      DATA TAU/.1D0 / , TAU0/1.D0/ , DEL0/1.D0/
      DATA (GUNIT(1,J),J=2,7)/6*1/,(GUNIT(3,I),I=2,7)/3*1,3*-1/,
     F     (GUNIT(2,I),I=2,7)/1,2,3,1,2,3/
     F     ,((GUNIT(I,J),I=1,3),J=0,1)/-1,0,0,-1,0,0/
      END
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX,F1,F2,F3,G1,G2,G3
      ICF=ICF+1
      F3=X(3)+.3D-1
      F2=F3+X(2)
      F1=F2+X(1)
      G1=F2+.9D-1* X(1)
      G2=F3+.7D-1* X(2)
      G3= .13D0 * X(3)+.3D-1
      FX=-(8.20437D3*DLOG(DABS(F1/G1))+ 9.00872D3*DLOG(DABS(F2/G2))
     1     +9.33046D3*DLOG(DABS(F3/G3)))
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),GRADF(NX),F1,F2,F3,G1,G2,G3
      ICGF=ICGF+1
      F3=X(3)+.3D-1
      F2=F3+X(2)
      F1=F2+X(1)
      G1=F2+.9D-1* X(1)
      G2=F3+.7D-1* X(2)
      G3= .13D0 * X(3)+.3D-1
      GRADF(1)=-8.20437D3* (1.D0/F1 - .9D-1/G1)
      GRADF(2)=-8.20437D3*(1.D0/F1-1.D0/G1)-9.00872D3*(1.D0/F2-.7D-1/G2)
      GRADF(3)=-8.20437D3*(1.D0/F1-1.D0/G1)-9.00872D3*(1.D0/F2-1.D0/G2)
     1         -9.33046D3*(1.D0/F3-.13D0/G3)
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      HXI= X(1)+X(2)+X(3)-1.D0
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      DO      100     J=1,3
      GRADHI(J)=1.D0
  100 CONTINUE
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GXI
      IF ( GUNIT(1,I+NH) .EQ. -1 ) CRES(I+NH)=CRES(I+NH)+1
      IF(I .LE. 3)      GXI=X(I)
      IF(I .GT. 3)      GXI=1.D0- X(I-3)
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX) ,GRADGI(NX)
      DO      100      J=1,3
      GRADGI(J)=0.D0
  100 CONTINUE
      IF(I .LE. 3)      GRADGI(I)=1.D0
      IF(I .GT. 3)      GRADGI(I-3)=-1.D0
      RETURN
      END


